Muon Spin Spectroscopy evidence of the violation of Anderson condition in magnetite. 
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We present new muon spectroscopy data on a Fe304 single crystal, revealing different spin pre- 
cession patterns in five distinct temperature ranges. A careful analysis of the local field and its 
straightforward modeling obtains surprisingly good agreement with experiments only if a very spe- 
cific model of localized charges violating Anderson condition, and a correlated muon local dynamics 
are implemented. Muon evidence for fluctuations just above the Verwey temperature, precursor of 
the low temperature charge localized state, is provided. 



The real nature of the charge ordered state in 
many transition metal oxide, notably manganites^ and 
magnetite^ ^, is still controversial. Do they correspond 
to distinct integer localized cation charges or to a much 
smaller charge disproportionation? The issue is of course 
also relevant to the nature of the carriers in the highly 
spin polarized metallic regimes of these oxides, which are 
of great interest for potential spintronic applications. 

The half-metallic character of magnetite, i.e. the fact 
that majority and minority-spin subbands are partially 
and completely filled, respectively, was re cently pointed 
out^ and attracted renewed interest"^"^"^"^ to this proto- 
typic magnetic material. In the original Verwey mode^, 
whose validity was not questioned until recentljE'^'^, it 
is due to charge delocalization between equal fractions of 
Fe^+ and Fe^+ ions, occupying the octahedral B sublat- 
tice of the inverse spinel structure, AB2O4, whereas Fe^+ 
ions stably occupy the tetrahedral A sublattice. The 
metal-insulator transition takes place at Ty (above 120 
K in good samples^^). The A and B spinel sublattices 
also correspond to the ferrimagnetic sublattices (Fig. [T]), 
hence all the B site spins are fully aligned, with easy axis 
along [111] above Ty. Below Ty the spin is parallel [001]. 
The B-sites form a pyrochlore network of corner-sharing 
tetrahedra (henceforth the B phyrochlore units, Bpcu). 

Assuming, with Verwey, localized Fe^+ and Fe^+ at 
B sites in the insulating state, AndersorP^ noticed that 
the rather low Ty implies a small activation energy. He 
proposed that instantaneous local charge configurations, 
both above and below Ty, must satisfy a condition of 
minimum local energy, where each Bcpu contains two 
Fe^+ and two Fe^+ ions. 

However both magnetic moments^ (Tab. [l| and local- 
ized charges below Ty do not correspond experimentally 
to those expected of Fe^+ and Fe^+. Bond valence sums^ 
and resonant x-ray scatterin^I^ indicate that the charge 
contrast among the B -site cations cannot be more than 
0.2e. Joint refinements^ of x-ray and neutron diffrac- 
tion data in the approximate P2/c symmetry relaxes the 
Anderson condition, supporting a model (hereafter, the 




FIG. 1: (Color online) spinel unit with A and B Fe ions and 
the T > Tv spin structure (left); low charge Bpcu (right) 
from Ref.[I3]with red O, blue Fe+2-^, green Fe+^ ^ 



Wright model) with Fe^"^+ and Fe^'^+ ions on the B sites. 
It corresponds to a [001] charge density wave (CDW, 
Fig. 7 in Ref. [T3|) composed of Bpcu with ¥e^-^"^/¥e^-^"^ 
occupancy ratio of 1:3 and tetrahedra with ratio 3:1 (the 
latter is shown in Fig. [T]). 

Direct evidence of a large symmetry reduction below 
Ty comes also from NMR^^ which identifies sixteen 
inequivalent B cations. Nuclear magnetic resonances also 
reveals a distinct spin reorientation transition^^, taking 
place slightly above Ty, at 126 K. More details on 

this aspect are given in Ref. 17, 

The reduced symmetry of the local environment 
may be probed directly by the muon spin precession 
around the magnetic field at the implantation site. We 
performed /iSR experiments on a high quality s ingle 
crystaP^at the Paul Scherrer Institut. Experiment^^SEIl 
with previous generation facilities and much lowe statis- 
tics measured only one precession frequency above Ty, 
and none below. Full muon site assignment and a rather 
detailed picture for charge localization emerge from our 
time dependent muon asymmetries in zero applied mag- 
netic field, which were already partially published with a 
very preliminary analysis^^. 

Fig. [2] shows a few representative precession pat- 
terns of the muon asymmetry, between 1.6 K and 300 
K. They are fitted to a sum of relaxing precessions. 
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FIG. 2: (Color online) Muon asymmetry in zero external field 
at selected temperatures, with best fit. 



A{t) = 2^oEj/jCos(27r75^^t)exp(-tV|/2)/3 where 
7 = 135.5 MHz/T is the muon magnetogyric ratio, Ao 
the total initial muon asymmetry, 5^^, crj/27r7, respec- 
tively, the local field intensity and its second moment for 
each transverse muon fraction fj. Additional non pre- 
cessing terms account for local field components parallel 
to the initial muon spin direction, which, in cubic and 
pseudocubic symmetries amount to an initial asymmetry 
of Ao/3. Further details and a color contour map of the 
Fast Fourier Transform (FFT) of the muon asymmetry 
are providecP^. 

Fig. [3] shows the temperature behavior of the local field 
strengths B^j^ where five distinct intervals are observed: 

i) For T>Tii = 250 K, S \\ [111], only one field, = 
0.42 T, is detected; the best fit has two components, one 
with fast (red diamonds) and one with slow relaxations 
(blue squares); 

a) For T/ = 160 K < T < T//, S \\ [111], two fields 
are detected, B^i ^ 0.36 T (green triangles) and B^2 ^ 
0.43 T (blue squares), with /2//1 = 3; red diamonds 
correspond again to the same field ^ 0.43 T, but 
with faster relaxation; 

in) For Tr < T < Tj a. third extra field B^s is detected, 
decreasing with temperature from roughly {B^i-\- B^i) /2 
towards an extrapolated value of 0.25 T, with fractions 
/i/(/2 + /s) = 3; 

iv) A sharp change takes place around Tr = 126(1) K, 
where, following the spin reorientation^^ from S \\ [111] 
(T > Tr) to S II [001] (T < Tr), B^i and B^2 collapse 
into B'^-^ = 0.435 T, while still survives; 

v) For T < Tv three fields are detected , B'^i^B'^2 ^ 
0.75 T and B'^^ ^ 1.06 T, with comparable fractions 
?^ 1/3 and large relaxation rates ai < (T2 < cr^] 

We concentrate here just on the field intensities, 
whereas finer details, such as relaxations, will be pub- 
lished elsewher^^. 
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FIG. 3: (Color online) Temperature dependence of the de- 
tected local muon fields, B^] the inset is a blow-up for 
T > Tv. 



TABLE I: Magnetic moments, in /xb, from Ref. [2] (cubic cell) 



Temperature 


FCA 


Fes 


S direction 


T>Tv 
T <Tv 


-4.20 
-4.44 


3.97 
4.17 


[111] 

[100] 



A unique muon interstitial site assignment justifies all 
these features. Its identification is first independently 
obtained from a simple electrostatic point-charge calcu- 
lation based on the notion that muons bind to oxygeii^^^ 
with bond lengths approximately equal to r^=l.l A. The 
minima of the electrostatic potential constrainecP^ 
on the sphere S{r^) centered on oxygen are shown in 
Fig.[4j3. Three equivalent minima (labeled a, 6 and c) are 
connected by a low potential path, separated by shallow 
barriers. They form a network in the lattice, as shown in 
Fig. 114. 

For an unmagnetized sample in zero external field the 
total local magnetic induction at a specific muon sit^^H 



(1) 



where is given by dipolar sums within a Lorentz 
sphere, with the known magnetic moments (Table |l|, 
Bl is the contribution from the Lorentz counter-sphere 
{Bl{T) = ^M(T), with M(0) equal to domain magne- 
tization, Bl(0) = 0.21 T), and Bhf an a-priori unknown 
isotropic^^ hyperfine contribution, also parallel to the do- 
main magnetization M. The three muon sites around 
each oxygen are crystallographically equivalent, hence 
they experience the same hyperfine field, but the elec- 
tron magnetic moment direction breaks the symmetrjB^^ 
yielding distinct dipolar fields. 

We can now discuss our experimental results, starting 
from high temperatures, T > T//, where one value of 
l^^l is detected. A single value implies^^ that the muon 
must be hopping among all equivalent sites. Since the 
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FIG. 5: (Color online) Sum of Bl + Buf (solid circles) with 
best fit to a power law (see text). 



FIG. 4: (Color online) 1) and 2) white balls are /x sites (a, 6, c) 
around a sphere S (color gradient) of radius 1.1 A, centered 
on oxygen (red online), on which (/)e(r) is calculated; Fca and 
Fes (yellow online). Spin orientation distinguishes 2), three 
times as frequent as 1). 3) Map of the electrostatic potential 
(j)e on <S, showing three equivalent minima (a, 6, c) with low 
valleys in between. 4) Muon site network seen from [Oil] - O 
ions lie in the (Oil) plane. 



spin orientation at the B site, S \\ [111], distinguishes two 
families of O ions, as shown in Fig. |4]l and 2, six large 
distinct local fields (j = 1,2 and a = a, 6, c) are predicted 
and no value of the hyperfine field Bhf can reconcile all 
of them with the experimenlF^ without invoking muon 
hopping. 

Unconstrained fast hopping above 250K takes place 
in other transition metal oxides, such as orthoferrite^^, 
cuprates^^ and manganite^. Its effect is to average out 
the dipolar fields, Ea=a,6,c[^d^ + 3Bj2]/4 = 0, thanks 



to cubic symmetry, yielding the same = Bl -\- Bhf 
for all muons, whence we obtain Bhf ^ 0.21 T. 

The same value of Bhf agrees with the experiment also 
for Tj < T < Tji^ if one assumes that now muon diffu- 
sion is restricted to fast tunneling among local a, b and 
c minima. This assumption, quite natural in view of 
the shallow barriers of Fig. [4]3, yields two average fields 
^/ij = Ea=a,6,c^M«j' j = 1,2, with moduH B^i = 0.36 
T, B^2 = 0.43 T and fractions in the ratio /2//1 = 3, as 
it is indeed observed in Fig. |3] (squares and triangles). 

The merging of Bj^i and into B'^^ below Tr, where 
the spin reorients, is also justified by the same assump- 
tions, since for || [100] (T < Tr) all oxygen ions become 
equivalent in the magnetic cell, yielding the same three 
dipolar field values in the three minima a, b and c, and 
their local average vanishes by cubic symmetry. Hence lo- 
cal tunneling predicts the average field B'^i = Bhf ^ Bl 
for all muons, in agreement with observation (squares 
in Fig. l3|. This same field value, B'^^ — Bhf -\- Bl^ 
is detected also below Ty, down to T=0. The same 
quantity may also be computed from {3B^2 + ^/ii)/4 for 



Tfi < T < Til and directly from B^ for T > T//. It is 
plotted versus temperature in Fig. |5) together with a fit 
to the power law B{T) = Bo(l -T/Tat)^. By imposing^^] 
Tat = 858 K, we obtain Bq = 0.447 T and /3=0.22 (the 
last parameter may be inaccurate, since T/T^ is lim- 
ited to 0.34). The overall agreement firmly establishes 
our site assignment and dipolar calculations, validating 
the simple electrostatic criterion and the two-stage muon 
diffusion. 

Let us consider now the range T < Ty. The ^e(^) 
minima may become inequivalent, depending on the local 
charge configuration (LOG) of the nearest Bpcu (Fig. [T]). 
Therefore, in agreement with higher temperature find- 
ings, we assume that muons either tunnel among equiv- 
alent minima or reside at inequivalent one. 

In the Verwey model, respecting Anderson condition, 
three distinct low charge muon LCC are identified, all of 
which with one lowest inequivalent (j)e minimum, hence 
no possible three-site tunneling. In the Wright model 
muons favour CDW charge troughs, located in specific 
[00/] planes. The oxygen ions in these planes correspond 
to those labeled Oi and O2 in Fig. [l] right; notice that 
Oi has three equivalent 0e minima, since its three nearest 
neighbor (n.n.) B ions are all Fe^*^. This leads to the 
correct prediction of a full local muon tunneling, hence 
of a low local field B'^^ = BhL even a T = 1.6 K. 

Table [IT] summarizes our findings for the two models, 
labeling each field value by the LCC of the B sites nearest 
neighbor to the muon. The agreement is very good with 
the Wright model, very poor with the Verwey model. 

The essential feature of a CDW along the c axis is the 
perfect correlation that it provides between the direction 
of the 0-Fe^"^ bond seen by the muon and that of the 
magnetic moments: they are all parallel to (001). This 
correlation grants the agreement with muon experiments 
and it is totally lost in the Verwey model. 

Finally, let us go back to the spin reorientation at Tr, 
around which a third, strongly temperature dependent 
field, 5^3, is observed. Its smooth decrease towards 0.2 
T for T ^ Ty must arise from fast fluctuations among 
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TABLE II: Total field intensity, in Tesla, at T = K for 
Verwey and Wright models (see text) for each LCC on B sites, 
labeled by three n.n. Fe valences, in bold. Overline indicates 
average among three sites, asterisk between two sites only. 



Model 


Local Charge Configurations 


Verwey 


223 


232 322 


Fields (T) 


1.14 0.70* 


1.85 0.98* 


Wright 


2.4 2.4 2.6 


2.4 2.4 2.4 


Fields (T) 


1.11 0.71* 


045 


Exp. Fields (T) 


1.06 0.80 


0.447 



two distinct local field configuration, with temperature 
dependent relative probabilities. We tentatively identify 
the two configurations as due to different local spin ori- 
entations. Indeed for a muon inside a [001] domain bub- 
ble the dipolar field vanishes. A static bubble within a 
larger [111] domain would also determine a cancellation 
of the Lorentz field, since a roughly equal, but oppo- 
site term is provided by the boundaries of the bubble 
itself. Hence the local field predicted by Eq. [l] would be 
B^3 = Bhf ^ 0.21 T. If the [001] bubble is fast fluctu- 
ating in a [111] background the muons may experience 
a temperature dependent average between the two static 
values, 0.21T and B^i (or 5^2)- 

This simple model, therefore, brings forward the fol- 
lowing picture: below T// = 150K there are regions 



where tiny bubbles of [001] spin orientation appear on 
a time-scale F"^ <C 20 ns, much shorter that the muon 
precession period. They are also characterized by a short 
coherence length, ^ (the bubble radius). If F"^ increases 
as Tv is approached, it shifts the weight in the muon 
average field B^s towards Bhf. This situation may well 
survive also below Tr, where a similar picture applies 
with exchanged roles: small bubbles of [111] inside a [001] 
domain also provide the cancellation of Bl. 

In conclusion we determine the muon location in Fe304 
and we detect a muon motion partially correlated with 
charge localization. Below Ty our findings are naturally 
reconciled with the structural model of Wright et al.^ pro- 
viding strong support for the violation of Anderson con- 
dition. We confirm a spin reorientation transition below 
126 K, from S ||[111] to S' || [001], precursor to the Verwey 
transition, and around it we deduce an inhomogeneous^, 
dynamic phase separation in fractions of the sample. 
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